function[u0]=pdeic(x)

u0=[sin(pi*x);sin(pi*x)];
%u0=[1;0];
%u0=[0;1];

